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There  are  mixed  opinions  on  the  suitability  of  dental  calculus  for  paleodietary  reconstruction  using  stable 
isotope  analysis.  We  examine  813C  and  515N  values  of  calculus  samples  from  two  regions,  central  Cali¬ 
fornia  in  the  USA  and  Sai  Island  in  the  Sudan.  When  atomic  C/N  ratios  are  less  than  12  in  calculus,  results 
show  positive  correlations  at  both  the  regional  and  individual  level  between  stable  isotopes  of  bone 
collagen  and  calculus,  suggesting  these  materials  track  similar  dietary  behaviors.  Correlations  are  still 
positive  but  lower  between  513C  values  of  calculus  and  bone  apatite.  Stable  isotope  ratios  of  calculus 
show  between  30%  and  50%  greater  variation  than  bone,  are  typically  enriched  in  15N  (mean  =  2.1  %o 
higher),  and  are  depleted  in  13C  relative  to  bone  collagen  (mean  =  0.8%o  lower)  and  apatite  (mean  =  6.4%o 
lower).  Calculus  from  multiple  teeth  was  analyzed  separately  for  seven  individuals  to  examine  intra¬ 
individual  variation.  Results  show  that  within  an  individual  S13C  varies  up  to  1.8%o,  and  815N  up  to 
2.1  %o,  which  may  explain  some  of  the  weak  bone-calculus  correlations  previously  reported  in  the 
literature.  When  atomic  C/N  ratios  are  greater  than  12,  calculus  correlates  more  poorly  with  bone 
collagen,  suggesting  these  samples  should  be  treated  with  caution. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Archaeological  teeth  often  contain  adhering  dental  calculus,  but 
this  material  is  often  overlooked  by  researchers  as  a  material  for 
archaeometric  analysis  (Lieverse,  1999).  Only  recently  has  the  po¬ 
tential  of  calculus  been  highlighted  in  research,  most  notably 
through  the  extraction  of  plant  starch  grains  from  the  calculus  of 
archaeological  humans,  and  more  recently  Neandertals  (Hardy 
et  al.,  2012;  Henry  et  al.,  2010).  Other  recent  analyses  focus  on 
recovering  bacterial  DNA  in  calculus  (Preus  et  al.,  2011),  charac¬ 
terizing  organic  compounds  by  chromatography  (Buckley  et  al., 
2014),  and  measuring  the  concentration  and  distribution  of 
particular  elements  such  as  calcium,  phosphorus,  and  manganese 
(Charlier  et  al.,  2010;  Klepinger  et  al.,  1977). 

Unlike  faunal  and  botanical  remains  from  archaeological  mid¬ 
dens,  which  inform  on  the  behaviors  of  groups  of  people,  the 
analysis  of  calculus  has  the  potential  to  inform  on  individual 
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behavior.  Such  data,  by  extension,  can  then  be  compared  between 
the  sexes,  between  individuals  of  different  statures,  and  compared 
to  health  and  social  status,  among  other  interesting  sub-groupings 
of  people  within  populations.  As  such,  there  exists  the  potential  to 
link  dietary  behaviors  to  a  wide  range  of  biological  and/or  cultural 
categories  in  ways  that  other  approaches  often  cannot.  Such 
individual-level  data  are  especially  relevant  in  evolutionary  an¬ 
thropology,  which  focuses  on  the  behaviors  of  individuals,  and  aids 
in  the  reconstruction  of  life  histories. 

Dental  calculus  is  composed  of  mineralized  plaque  (Hillson, 
1996).  It  consists  of  a  range  of  inorganic  minerals,  including  hy¬ 
droxyapatite,  octacalcium  phosphate,  and  whitlockite,  but  can  also 
encompass  a  suite  of  organic  components,  including  plant  fibers, 
glycoproteins,  proteins,  lipids,  and  carbohydrates  (Kani  et  al.,  1983; 
Lieverse,  1999).  Deposits  of  silicates,  including  quartz,  calcium  sil¬ 
icates,  and  calcium-aluminum  silicates  can  also  precipitate  into  the 
porous  structure  of  calculus  in  archaeological  samples  when  teeth 
are  in  contact  with  ground  water  (Klepinger  et  al.,  1977).  In  clinical 
studies,  supragingival  calculus  most  commonly  forms  on  the 
lingual  surfaces  of  anterior  mandibular  teeth  and  the  buccal  sur¬ 
faces  of  maxillary  molars  adjacent  to  salival  ducts  (White, 
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1997:513).  Progression  can  begin  after  tooth  eruption,  and  may  take 
a  decade  or  more  to  form.  In  teeth  that  are  not  regularly  cleaned  of 
plaque,  calculus  can  peak  in  volume  around  30  years  of  age  (White, 
1997:510).  The  organic  component  provides  the  source  material  for 
stable  carbon  and  nitrogen  isotope  analysis  of  dental  calculus 
(Poulson  et  al,  2013). 

Stable  isotope  analysis  of  collagen  and  apatite  from  bone  and 
teeth  provide  important  windows  on  our  understanding  of  ancient 
diet  and  mobility  practices.  Such  analyses  have  highlighted  differ¬ 
ences  and/or  similarities  in  diet  of  elites  vs.  commoners  in  complex 
societies  (e.g.,  Ambrose  et  al.,  2003;  Somerville  et  al.,  2013; 
Voutsaki  et  al.,  2013),  differential  mobility  practices  among  adults 
suggesting  a  societal  preference  for  matrilocal  or  patrilocal  post- 
marital  residence  (e.g.,  Dupras  and  Schwarcz,  2001;  Scheeres 
et  al.,  2013;  Toyne  et  al.,  2014),  and  age  at  weaning  and  dietary 
differences  by  sex  (e.g.,  Dupras  et  al.,  2001;  Dupras  and  Tocheri, 
2007;  Eerkens  and  Bartelink,  2013). 

A  recent  analysis  by  Scott  and  Poulson  (2012)  employed  stable 
isotope  analysis  of  cleaned  but  otherwise  untreated  bulk  dental 
calculus  samples.  They  showed  that  stable  isotope  signatures  in  the 
calculus  were  within  the  range  of  what  would  be  expected  of  bone 
collagen,  though  they  did  not  analyze  bone  collagen  from  the  same 
set  of  individuals.  From  this,  they  concluded  that  calculus  may  be  a 
suitable  material  for  reconstructing  paleodiet.  In  a  second  study, 
Poulson  et  al.  (2013)  examined  stable  carbon  and  nitrogen  isotopes 
in  ancient  dental  calculus  from  northern  Chile,  and  found,  again, 
that  results  matched  what  was  expected  from  bone  collagen.  For 
seven  of  the  28  burials  studied,  they  examined  calculus  from  two 
different  areas  of  the  dentition  to  evaluate  intra-individual  varia¬ 
tion  (Poulson  et  al.,  2013:4578).  However,  as  in  Scott  and  Poulson 
(2012),  their  sample  did  not  include  isotopic  data  on  paired  calcu¬ 
lus  and  collagen  from  the  same  individuals. 

Contrasting  these  results,  Salazar-Garcia  et  al.  (2014)  concluded 
that  513C  vs  815N  values  of  dental  calculus  showed  little  correlation 
to  values  from  bone  collagen  from  the  same  individuals  in  a  sample 
of  35  burials  from  Medieval  Spain.  While  calculus  and  bone  813C 
and  515N  values  were  found  to  be  similar  for  some  individuals,  for 
others  they  differed  up  to  10%o  in  813C  and  3.5%o  in  815N.  As  a  result, 
they  did  not  recommend  stable  isotope  analysis  of  calculus  as  a 
means  to  reconstruct  paleodiet. 

Given  these  contrasting  results,  we  wanted  to  re-evaluate  the 
question  of  using  calculus  stable  isotope  analysis  as  a  proxy  for  diet 
by  conducting  additional  analyses  on  paired  calculus  and  bone 
collagen  samples,  in  this  case  from  western  North  American 
hunter-gatherers  and  northeastern  African  agriculturalists.  We  also 
examine  correlations  with  813C  in  bone  apatite  for  a  subsample  of 
these  individuals,  on  which  no  previous  research  has  been  con¬ 
ducted.  Finally,  like  Poulson  et  al.  (2013),  we  extracted  and  analyzed 
separately  calculus  samples  from  multiple  teeth  for  seven  in¬ 
dividuals,  to  estimate  the  degree  of  intra-individual  variation, 
which  we  compare  to  intra-site  and  inter-site  variation. 

2.  Samples 

Supragingival  calculus  samples  were  taken  from  46  individuals 
from  two  disparate  regions.  All  sites  were  occupied  during  the 
latter  half  of  the  Holocene.  Table  1  shows  the  samples  included  in 
the  analysis,  listed  by  site. 

The  North  American  sample  includes  36  individuals  from  five 
Middle  and  Late  Holocene  sites  in  central  California  (see  Fig.  1). 
Sites  are  included  from  a  range  of  environmental  zones,  from  the 
southern  San  Francisco  Bay  estuary,  to  brackish-water  Suisun 
Marsh,  to  more  inland  locations,  in  an  attempt  to  capture  a  range  of 
potential  isotope  values.  Although  not  all  the  burials  included  in 
this  analysis  were  directly  dated,  all  the  cemeteries  have  multiple 


Table  1 

Distribution  of  samples  included  in  this  study. 


Geographic 

region 

Site 

Calculus 

only 

#  calculus- 

collagen 

pairs 

#  calculus- 

apatite 

pairs 

#  indiv.  with 

multiple 

calculus 

California,  USA 

CA-ALA-554 

0 

4 

0 

0 

California,  USA 

CA-CCO-548 

1 

6 

6 

2 

California,  USA 

CA-SOL-270 

0 

21 

8 

0 

California,  USA 

CA-SOL-11 

0 

2 

2 

0 

California,  USA 

CA-SCL-134 

0 

1 

1 

0 

California,  USA 

CA-SCL-287 

0 

1 

1 

0 

Republic  of  the 
Sudan 

Sai  Island 

5 

5 

5 

5 

Total 

6 

40 

23 

7 

bone  collagen  radiocarbon  dates  indicating  that  each  was  in  use 
during  a  relatively  narrow  window  of  time  (see  Eerkens  et  al.,  2013; 
Estes  et  al.,  2012;  Leventhal  et  al.,  2010).  The  sites  (with  number  of 
burials  and  age  range)  include  CA-ALA-554  (n  =  4; 
2000-600  cal  BP),  CA-CCO-548  (n  =  7;  4000-3000  cal  BP),  CA-SOL- 
270  (n  =  21;  2500-1700  cal  BP),  CA-SOL-11  (n  =  2; 
2300-1700  cal  BP),  CA-SCL-134  (n  =  1;  2800-2300  cal  BP),  and 
CA-SCL-287  (n  =  1 ;  2000-1300  BP).  Previous  analyses  indicate  the 
individuals  represented  in  this  study  were  hunter-fisher-gatherers, 
relying  on  a  wide  range  of  plant  and  animal  resources.  Data  indicate 
that  all  cemeteries  are  associated  with  relatively  permanent  (i.e., 
sedentary)  villages.  For  two  individuals  from  CCO-548,  multiple  cal¬ 
culus  samples  were  taken  from  different  teeth.  In  total,  39  calculus 
samples  were  analyzed  from  California. 

The  African  sample  includes  ten  individuals  from  a  single  site  on 
Sai  Island  in  the  Republic  of  the  Sudan.  Sai  Island  is  the  largest  is¬ 
land  of  the  Nile  valley  located  between  the  2nd  and  3rd  cataracts  of 
Nubia  (Fig.  2).  The  ten  individuals  were  buried  in  an  elite  necropolis 
belonging  to  the  Meroitic  Kingdom  (2400  BP  to  1600  BP;  de  Voogt 
and  Francigny,  2012),  though  the  skeletal  material  for  this  study 
dates  to  a  more  narrow  temporal  window  between  the  first  and  the 
third  century  CE  (1950-1750  BP).  The  Meroites  from  the  northern 
part  of  the  Nubian  kingdom,  where  Sai  Island  is  located,  mostly 
relied  on  riverine  agriculture  supported  by  the  alluvium  deposited 
during  the  annual  flood.  Herding  seems  to  have  been  widespread 
but  very  limited  compared  to  the  southern  half  of  the  kingdom  and 
the  capital,  where  rainland  agriculture  was  also  developed.  Though 
archaeobotanical  data  are  limited,  it  is  generally  assumed  that 
millet,  sorghum,  wheat  and  barley  were  cultivated  by  the  Meroites 
(Thompson  et  al.,  2008;  Welsby,  1998:155-161).  Bone  samples 
were  only  available  for  five  individuals.  In  addition,  for  five  in¬ 
dividuals,  multiple  calculus  samples  from  different  teeth  were 
analyzed.  In  total,  26  calculus  samples  were  analyzed  from  Sai 
Island. 

3.  Methods 

Calculus  preparation  methods  followed  Scott  and  Poulson 
(2012).  Samples  were  removed  from  the  tooth  surface  using  a 
scalpel,  and  were  cleaned  but  not  chemically  treated.  Calculus  was 
rinsed  in  deionized  water  and  non-calculus  material,  such  as 
sediment,  was  removed  by  hand  with  a  set  of  tweezers  under  a 
low-power  microscope.  Calculus  was  then  powdered  in  a  mortar 
and  pestle.  For  individuals  with  multiple  calculus  samples,  the 
resulting  813C  and  815N  values  were  averaged  (separately  for  C  and 
N)  to  compare  to  bone  collagen  and  apatite  values. 

Between  1  and  12  mg  of  powdered  calculus  was  submitted  to 
the  Stable  Isotope  Facility  at  UC  Davis,  and  measured  for  total  C  and 
N  (in  pg),  and  813C  and  815N,  using  a  PDZ  Europa  ANCA-GSL 
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Fig.  1.  Map  of  central  California  showing  locations  of  sites  included  in  this  study. 


elemental  analyzer  interfaced  to  a  PDZ  Europa  20-20  isotope  ratio 
mass  spectrometer  (Sercon  Ltd.,  Cheshire,  UK).  For  samples  with 
sufficient  C  (ca.  100  pg),  the  long-term  standard  deviation  on  this 
instrument  is  0.2%o  for  513C,  and  for  samples  with  greater  than 
20  pg  N,  the  long-term  standard  deviation  is  0.3%o  for  815N.  In¬ 
strument  precision  decreases  for  smaller  samples,  and  we  only 
report  values  for  samples  containing  greater  than  20  fig  C  or 
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Fig.  2.  Map  of  northeast  Africa,  showing  location  of  Sai  Island. 


10  fig  N.  Analyses  were  limited  by  the  amount  of  calculus  present 
on  a  tooth  and  some  of  the  samples  in  this  analysis  weighed  less 
than  1  mg  but  were  submitted  nevertheless.  Several  of  these  pro¬ 
duced  enough  C  but  not  enough  N  for  analysis.  As  well,  for  some 
individuals  in  the  California  sample,  we  combined  calculus  from 
more  than  one  tooth  to  generate  enough  material  for  stable  isotope 
analysis  (see  Supplementary  information). 

Bone  collagen  extraction  for  all  samples  followed  a  modified 
Longin  procedure  (Longin  1971).  Approximately  1-2  g  of  cortical 
bone  was  cleaned  of  any  surface  contamination  by  drilling  exposed 
surfaces  with  a  diamond-studded  bit  and  then  sonicating  the 
sample  in  deionized  H2O.  The  sample  was  demineralized  with  a 
solution  of  .5  M  hydrochloric  acid  (HC1).  The  bone  was  then  washed 
and  soaked  in  .125  M  NaOH  (sodium  hydroxide)  to  remove  humic 
acids.  The  sample  was  then  placed  in  slightly  acidic  water  in  a 
70-90  °C  oven  for  approximately  24  h  to  solubilize  collagen  and 
freeze  dried  to  remove  the  water  and  isolate  the  collagen  fraction. 
While  the  Californian  collagen  samples  were  analyzed  on  the  same 
mass  spectrometer  as  the  calculus  samples  above,  the  Nubian 
collagen  samples  were  analyzed  on  a  Delta  V  Advantage  Mass 
Spectrometer  with  a  Carlo  Erba  NC2100  Elemental  Analyzer  at  the 
Colorado  Plateau  Stable  Isotope  Laboratory  at  Northern  Arizona 
University. 

Bone  apatite  preparation  started  with  the  same  cleaned  bone, 
which  was  powdered  in  an  agate  mortar  and  pestle.  Organics  were 
removed  by  adding  a  1.5%  sodium  hypochlorite  at  a  ratio  of  0.04  ml 
solution/mg  sample  (Koch  et  al.,  1997),  followed  by  a  diagenetic 
wash  composed  of  a  1  M  acetic  acid  solution  (at  the  same  ratio  of 
.04  ml  solution/mg  sample).  The  sample  was  then  rinsed  with  dE^O 
and  dried.  Apatite  samples  from  California  were  analyzed  on  a  GVI 
Optima  Stable  Isotope  Ratio  Mass  Spectrometer  at  the  Stable 
Isotope  Lab  in  the  Geology  Department  at  University  of  California, 
Davis,  with  513C  measured  on  the  carbonate  component  of  bone 
apatite.  Bone  apatite  for  the  Nubian  samples  were  prepared  using 
the  same  method  as  previously  described,  and  analyzed  on  a 
Thermo-Finnigan  Gasbench  II  and  a  Thermo-Finnigan  Delta  Plus 
Isotope  Ratio  Mass  Spectrometer  at  the  Colorado  Plateau  Stable 
Isotope  Laboratory  at  Northern  Arizona  University. 
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4.  Results 

Complete  data  are  available  as  Supplementary  Materials  to  this 
article.  While  there  was  enough  C  (>20  pg)  in  all  calculus  samples, 
ten  samples  did  not  produce  enough  N  (>10  pg)  to  be  included  in 
the  analysis.  In  general,  there  was  a  strong  and  linear  correlation 
between  the  sample  weight  and  the  recovered  total  C  (r2  =  0.86) 
and  N  (r2  =  0.77),  suggesting  roughly  similar  elemental  composi¬ 
tion  across  the  majority  of  samples.  Overall,  calculus  samples 
averaged  50  pg  of  C  per  mg  of  calculus  (i.e.,  5%  of  calculus  by 
weight),  but  just  7  pg  of  N  per  mg  of  calculus  (0.7%  of  calculus  by 
weight).  These  figures  are  very  similar  to  Salazar-Garcia  et  al. 
(2014),  who  report  figures  of  4.7%  and  0.7%  (C  and  N  respec¬ 
tively),  and  Scott  and  Poulson  (2012)  at  4.9%  and  0.7%  for  Medieval 
Spanish  burials,  but  about  half  that  reported  by  Poulson  et  al.  (2013) 
for  Late  Holocene  burials  from  Chile  (11.7%  C  and  1.4%  N).  We  note 
that  there  was  some  variation  in  our  sample  by  region.  The  Nubian 
samples,  on  average,  produced  slightly  more  carbon  (mean  =  5.9%) 
than  the  California  samples  (mean  =  4.5%),  but  nearly  double  the 
nitrogen  (mean  =  0.94%  vs.  0.47%). 

Fig.  3  shows  the  correlations  between  isotopic  values  from  cal¬ 
culus  and  bone  with  symbols  used  to  denote  different  sites.  The  top 
and  bottom  panels  show  813C  and  815N  values  for  calculus  vs.  bone 
collagen,  respectively,  while  the  middle  panel  shows  513C  values  for 
calculus  vs.  bone  apatite.  Also  shown  are  best-fit  regression  lines 
and  Pearson's  r2  correlation  statistics  for  all  individuals  in  the  study 
together  (heavy  dashed  line),  as  well  as  regression  lines  for  indi¬ 
vidual  sites  (light  dotted  lines). 

The  data  points  in  Fig.  3  clearly  separate  the  Californian  from 
Nubian  samples,  for  both  calculus  and  bone.  Our  main  goal  here  is 
not  to  reconstruct  regional  or  individual  differences  in  dietary  be¬ 
haviors  from  calculus,  but  rather,  to  explore  patterns  in  calculus 
isotopes  relative  to  those  recorded  in  bone.  Nevertheless,  if  dental 
calculus  stable  isotope  signatures  reflect  diet,  as  bone  does,  these 
differences  suggest  quite  different  diets  in  the  two  continents.  Of 
course,  it  is  not  surprising  to  find  such  dietary  differences  in  sites 
representing  hunter-gatherers  in  California,  on  the  one  hand,  and 
agriculturalists  in  Nubia,  on  the  other.  However,  if  the  differences 
are  due  to  other,  non-dietary  factors  such  as  environment,  popu¬ 
lation  genetic  history,  or  the  composition  of  oral  bacteria  commu¬ 
nities,  we  would  expect  similar  differences  in  such  widely 
divergent  samples.  This  study  is  currently  unable  to  address  this 
issue,  but  note  that  the  different  sites  within  California,  populations 
with  more  similar  genetic  and  environmental  histories,  also  tend  to 
cluster  in  the  graph.  In  any  case,  these  data  point  to  unique  calculus 
isotopic  values  from  region  to  region  and  site  to  site,  and  that  these 
patterns  tend  to  mimic  those  in  bone  collagen.  We  note  three 
additional  patterns  below. 

First,  when  all  individuals  are  examined  collectively,  the  calcu¬ 
lus  isotopes  are  positively  correlated  with  bone  collagen  isotope 
values  for  both  513C  and  815N,  but  are  more  weakly  correlated  with 
bone  apatite  513C  values.  The  correlation  is  strongest  for  815N  (Fig.  3, 
lower  panel),  where  there  is  an  obvious  and  positive  relationship 
(r2  =  0.62).  The  correlation  between  calculus  and  bone  collagen 
813C  is  weaker  (r2  =  0.13),  due  primarily  to  elevated  813C  calculus 
values  from  one  site,  CCO-548.  If  the  CCO-548  samples  are 
excluded,  the  collective  correlation  for  813C  is  much  stronger 
(r2  =  0.61)  and  nearly  equal  to  the  correlation  with  815N.  On  the 
other  hand  the  figure  suggests  a  positive  but  weak  correlation  be¬ 
tween  calculus  and  bone  apatite  813C  (r2  =  0.03).  Removing  the 
CCO-548  samples  again  improves  the  correlation  with  bone  apatite 
813C  (r2  =  0.36),  but  it  is  still  notably  lower  than  correlations  with 
bone  collagen. 

Note  further  that  these  relationships  also  tend  to  hold  on  a  site- 
by-site  basis,  fable  2  shows  correlation  statistics,  and  differences  in 
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Fig.  3.  Comparison  of  513C  in  bone  collagen  (top)  and  apatite  (middle)  and  515N  from 
bone  collagen  (bottom)  against  isotopes  from  calculus,  by  individual. 


the  statistical  range  of  isotopic  values  for  calculus  vs.  bone  collagen, 
and  apatite,  for  sites  with  more  than  three  paired  measurements. 
The  strength  of  the  correlation  within  sites  is  lower,  as  measured  by 
r2,  but  still  positive  for  calculus  and  bone  within  each  site.  The  only 
exception  is  the  calculus-bone  apatite  813C  values  from  Sai  Island, 
where  there  is  no  correlation  (r2  =  0). 

Second,  calculus  isotopes  typically  show  greater  variation,  as 
measured  by  the  range,  than  bone  collagen.  This,  again,  holds 


Table  2 

Correlation  statistics,  and  the  magnitude  of  the  increase  in  range  for  calculus. 


Calculus  VS. 

SOL-270 

CCO-548 

ALA-554 

Sai  Island 

Overall 

Collagen 

r2 

0.01 

0.06 

0.23 

0.04 

0.13 

513C 

Range 

+2.4%o 

+4.8%o 

+0.2%o 

+0.6%o 

+2.8%o 

Apatite 

r2 

0.03 

0.56 

n/a 

0.00 

0.03 

S13C 

Range 

+0.4%o 

+2.2%o 

n/a 

+2.0%o 

+  1.1%0 

Collagen 

r2 

0.17 

n/a 

0.36 

0.02 

0.62 

515N 

Range 

+0.3%o 

— 0.2%o 

— 0.8%o 

— 0.8%o 

+2.4%o 

Notes:  n/a  =  no  data  available  or  data  set  too  small  for  correlation  (n  <  3). 
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within  sites  and  overall  (Table  2).  When  considered  as  a  collective 
data  set,  calculus  513C  has  a  2.8%o  greater  range  (46%  higher)  than 
bone  collagen  513C  and  1.1  %o  greater  range  (16%  higher)  than  bone 
apatite  513C.  Likewise,  calculus  815N  has  a  2.4%o  greater  range  (30% 
greater)  than  bone  collagen  815N.  Similarly,  when  we  examine  data 
within  sites,  calculus  513C  is  between  27%  (Sai  Island)  and  364% 
(CCO-548)  more  variable  than  bone  collagen  813C,  and  between  10% 
(SOL-270)  and  48%  (CCO-548)  more  variable  than  bone  apatite  513C. 
On  the  other  hand,  calculus  815N  is  less  variable  than  bone  collagen 
at  CCO-548,  ALA-554  and  Sai  Island,  but  is  more  variable  at  SOL- 
270. 

Third,  the  data  also  show  that  the  absolute  813C  values  in  cal¬ 
culus  tend  to  be  lower  (i.e.,  more  negative),  but  815N  elevated, 
relative  to  bone.  Thus,  individuals  tend  to  fall  in  the  lower  right  (i.e., 
below  the  diagonal)  in  the  plots  for  813C,  but  in  the  upper  left 
(above  the  diagonal)  for  815N  in  Fig.  3.  The  only  exception  to  this 
pattern  are  the  individuals  from  CCO-548,  where  calculus  813C  is 
elevated  relative  to  bone  collagen. 

Because  human  bone  collagen  has  a  relatively  fixed  chemical 
composition,  the  atomic  C/N  ratio  is  a  useful  indicator  of  sample 
quality  in  ancient  bone  collagen  (DeNiro,  1985;  Weber  et  al.,  2005). 
For  bone  collagen,  studies  show  that  samples  with  C/N  between  2.7 
and  3.3  produce  813C  and  815N  values  within  expected  ranges,  while 
C/N  falling  above  or  below  this  range  more  often  produces  unusual 
or  outlying  isotopic  values.  All  bone  collagen  samples  in  this  study 
produced  C/N  within  this  range.  No  such  accepted  C/N  measure 
exists  for  calculus,  and  the  more  variable  biomolecular  composition 
of  calculus  makes  it  difficult  to  predict  an  acceptable  range  for 
“good”  samples.  Calculus  samples  in  this  study  vary  from  29.0  to  5.9 
in  C/N,  with  a  mean  of  9.9  ±  3.2  (1  SD). 

However,  our  data  show  that  there  are  notable  patterns  in  C/N 
ratios  and  the  calculus-bone  offset.  Fig.  4  plots  C/N  in  calculus 
against  the  813C  value  of  calculus  minus  the  513C  value  of  bone 
collagen  (the  calculus-bone  collagen  offset).  For  comparison,  we 
also  include  data  published  by  Salazar-Garria  et  al.  (2014)  from  a 
medieval  period  site  in  Spain.  The  figure  shows,  first,  that  513C  is 
generally  lower  in  calculus  than  in  bone  collagen,  as  discussed 
above,  resulting  in  negative  values  on  the  y-axis.  Second,  the 
calculus-bone  collagen  513C  offset  is  increasingly  variable  for 
samples  where  the  calculus  C/N  ratio  is  greater  than  12.  One  site  in 
particular  (CCO-548)  has  consistently  high  C/N,  and  unlike  other 
sites,  the  calculus-bone  collagen  offset  for  813C  is  notably  positive 
for  nearly  all  samples  at  this  site  (see  Fig.  3,  top  panel).  For  calculus 
samples  with  C/N  less  than  12,  the  mean  calculus-bone  collagen 
offset  for  813C  was  -1.3%o  (0.9  =  1  SD).  Likewise,  for  these  samples 
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Fig.  4.  Calculus  -  Bone  collagen  offset  for  513C  vs.  calculus  C/N  ratio. 


the  mean  calculus-bone  apatite  813C  offset  was  -7.8%o  (1.6  =  1  SD) 
and  the  mean  calculus-bone  offset  for  815N  was  2.1  %o  (1.8  =  1  SD). 
Note  that  this  pattern  among  calculus  samples  with  high  C/N  ratios 
also  holds  among  the  Medieval  Spain  samples  reported  by  Salazar- 
Garcia  et  al.  (2014)  from  El  Raval,  and  may  explain  some  of  their 
reservations  about  using  calculus  to  predict  bone  collagen  values. 

Fig.  5  shows  the  813C  and  815N  results  of  repeated  sampling  of 
calculus  from  different  teeth  within  the  same  individual  for  two 
CCO-548  and  five  Sai  Island  burials.  The  x-axis  denotes  different 
individuals,  reproduced  in  the  same  order  for  the  upper  and  lower 
panes,  and  the  small  black  circles  represent  all  samples  from  that 
particular  site.  The  figure  shows  that  calculus  samples  from 
different  teeth  within  an  individual  do  not  produce  identical 
isotope  values,  but  differ  by  up  to  1.8%o  for  513C  and  2.1  %o  for  815N. 
This  is  slightly  greater  than  the  maximum  intra-individual  differ¬ 
ences  of  1.0%o  and  2.0%o  found  by  Poulson  and  colleagues  for  cal¬ 
culus  samples  in  Chile  (Poulson  et  al.,  2013).  As  well,  it  is  much 
greater  than  the  degree  of  instrument  precision  (0.2-0.3%o). 

At  the  same  time,  the  figure  also  shows  that  the  ranges  and 
standard  deviations  of  values  within  an  individual  are  less  than  the 
range  and  standard  deviation  of  values  within  each  site  as  a  whole. 
Thus,  the  standard  deviation  of  513C  values  in  all  Sai  Island  samples 
is  1.3%o,  while  standard  deviations  for  T035-individual  2  and  T030- 
individual  5  are  0.6%o  and  0.3%o,  respectively  (the  only  individuals 
with  more  than  three  samples).  In  spite  of  small  sample  sizes  and  a 
lack  of  information  about  the  underlying  statistical  distribution  of 
our  isotope  data,  two-sample  F-tests  comparing  equality  of  vari¬ 
ance  for  intra-tooth  vs.  intra-population  values  for  813C  are  signif¬ 
icant  for  both  individuals  (p  =  0.01  and  0.04,  respectively). 
Likewise,  the  standard  deviation  in  513C  for  all  CCO-548  samples  is 
7.2%o,  while  burial  59  has  a  standard  deviation  of  just  1.2%o  (F-test 
comparing  equality  of  variance  is  not  significant,  with  p  =  0.10,  but 
the  sample  size  is  much  smaller).  The  same  general  pattern  holds 
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more  weakly  for  S15N  values,  though  there  are  fewer  samples 
available  for  analysis,  precluding  performance  of  F-tests.  Thus, 
intra-individual  variation  for  515N  is  closer  to  the  inter-individual 
variation  for  a  particular  site.  While  we  do  not  place  great  inter- 
pretational  weight  on  the  results  from  the  statistical  tests,  they  do 
support  our  general  argument  that  the  degree  of  isotopic  variation 
is  less  within  individuals  than  within  the  population  as  a  whole. 

We  also  examined  isotopic  variation  by  tooth  type  and  position 
in  the  mouth.  All  but  one  of  our  repeat  samples  are  from  mandib¬ 
ular  teeth,  thus  we  are  unable  to  compare  teeth  from  the  upper  and 
lower  oral  cavity.  Teeth  from  the  left  and  right  side  of  the 
mandibular  dentition  overlap  completely,  suggesting  no  difference. 
However,  the  data  do  indicate  differences  between  posterior  and 
anterior  teeth.  For  example,  the  incisors  and  canines  of  T035-individual 
2  at  Sai  Island  have  slightly  more  negative  513C  values  (mean  =  -  18.5%o; 
n  =  4),  relative  to  premolars  and  molars  (mean  =  -17.6%o;  n  =  6).  Even 
though  the  sample  size  is  small,  a  t- test  comparing  calculus  from 
anterior  vs.  posterior  teeth  (2-tailed,  assuming  unequal  variance)  is 
statistically  significant  (p  =  0.02)  for  this  individual.  On  the  other 
hand,  three  other  individuals  show  the  opposite  pattern,  though 
the  differences  are  not  significant.  813C  values  of  calculus  from  a 
molar  from  T312  is  more  negative  (-21.3%o)  relative  to  the  canine 
(— 19.5%o),  calculus  from  the  molar  of  T034A-Individual  1  is  more 
negative  (-17.9%o)  relative  to  the  canine  (-17.4%o),  and  the  molar 
and  premolar  are  more  negative  (mean  =  -16.2%o)  than  the  incisor 
and  canine  (mean  =  -15.8%o)  from  T030-Individual  5.  There  are  not 
enough  815N  values  to  make  a  similar  comparison  for  nitrogen. 
These  results,  particularly  as  they  pertain  to  intra-individual  vari¬ 
ation,  are  discussed  further  below. 

5.  Discussion 

Our  results  suggest  several  items  of  interest  for  isotopic  analyses 
of  human  dental  calculus.  First,  calculus  C  and  N  isotope  values 
appear  to  conserve  the  between-site  differences  that  are  present  in 
bone  collagen  and  bone  apatite  isotopes.  This  is  consistent  with  the 
studies  of  Scott  and  Poulson  (2012)  and  Poulson  et  al.  (2013),  who 
found  that  stable  isotopes  of  calculus  were  similar  to  previously 
published  studies  of  bone  collagen.  In  support  of  those  studies,  we 
agree  that  isotopic  analysis  of  calculus  may  be  a  useful  tool  for 
estimating  general  dietary  differences  between  regions  or  over 
time  in  a  single  area  when  bone  is  not  available. 

Second,  while  preserving  inter-regional  signals,  there  appears  to 
be  an  offset  in  our  samples  between  calculus  and  bone  in  that  the 
former  are  often  depleted  in  13C  but  enriched  in  15N.  If  calculus 
isotopes  reflect  diet,  this  finding  is  consistent  with  the  notion  that  a 
significant  portion  of  the  carbon  in  calculus  is  derived  from  organic 
plant  matter,  such  as  starch  and  cellulose.  These  biomolecules  are 
high  in  carbon  but  low  in  nitrogen,  and  would  therefore  contribute 
mainly  to  the  813C  signatures  in  calculus.  Terrestrial  C3  plants 
typically  have  lower  513C  values  than  tissues  derived  from  fish  or 
other  animal  sources,  especially  proteins  which  contribute  the 
majority  of  carbon  to  bone  collagen  (Fernandes  et  al.,  2012).  This 
could  account  for  the  lower  813C  values  in  calculus  relative  to  bone. 

The  higher  mean  815N  values  in  calculus  (averaging  2.1  %o  higher 
in  this  study)  could  likewise  be  due  to  the  presence  of  significant 
amounts  of  bacteria  and  bacterial  byproducts  preserved  in  calculus. 
The  microbial  mass  preserved  in  calculus  is  expected  to  be  enriched 
in  15N  relative  to  the  proteins  ingested  by  humans  (and  routed  to 
synthesize  bone  collagen).  Indeed,  recent  studies  of  human  dental 
plaque  show  tremendous  bacterial  diversity,  with  well  over  300 
different  species  present  in  a  single  sample  (Paster  et  al.,  2001 ).  We 
are  not  aware  of  isotopic  studies  specifically  on  human  oral  bac¬ 
teria,  but  studies  on  other  bacterial  colonies  indicate  extremely 
complex  fractionation  processes  that  vary  greatly  depending  on  the 


energy  sources  introduced  (e.g.,  protein,  glucose,  cellulose,  etc.)  and 
interactions  between  different  bacterial  species  (Wattiaux  and 
Reed,  1995).  Variation  among  different  human  populations,  and/ 
or  individuals  within  populations,  in  the  species  diversity  and 
number  of  bacteria  in  the  oral  cavity  are  likely  to  contribute  to 
offsets  between  dietary  815N,  bone  collagen  815N,  and  calculus  815N. 

Our  results  also  show  that  the  magnitude  of  the  offset  in  cal¬ 
culus  relative  to  bone  can  vary  by  site  and  region.  Even  within 
California,  CCO-548  appears  to  have  elevated  813C  values  in  calculus 
relative  to  bone  while  other  sites  show  the  opposite  pattern.  This 
suggests  that  it  may  not  be  possible  to  simply  add  or  subtract  a 
given  amount  from  calculus  813C  and  815N  to  approximate  what 
would  be  observed  in  bone.  These  inter-site  differences  in  calculus- 
bone  offsets  point  to  either  variation  in  the  composition  of  calculus 
from  site  to  site,  which  may  be  indicative  of  differences  in  paleodiet, 
or  to  preservational  conditions  at  different  sites  in  which  case 
calculus  isotopes  are  less  useful  for  reconstructing  paleodiet. 
Lacking  additional  information  on  the  composition  of  calculus 
samples  and/or  a  good  index  for  estimating  the  state  of  preserva¬ 
tion  (as  exists  for  bone  collagen  using  collagen  yield  and  C/N  ratios), 
it  is  difficult  to  evaluate  which  of  these  scenarios  accounts  for  the 
patterns  reported  (though  see  below  for  discussion  on  C/N  ratios  in 
calculus). 

Third,  while  regional  signals  are  preserved  within  sites,  calculus 
isotope  values  only  marginally  preserve  the  inter-individual  vari¬ 
ation  in  C  and  N  isotopes  observed  in  bone  collagen  and  apatite. 
Correlations  between  bone  and  calculus  isotopes  are  positive 
within  sites,  suggesting  a  relationship  exists  between  the  two 
indices,  but  low,  indicating  that  the  nature  of  that  relationship  is 
more  complex.  Of  course,  this  finding  does  not  necessarily  mean 
that  calculus  isotopes  do  not  inform  on  inter-individual  differences 
in  diet.  As  we  have  shown,  there  is  significant  inter-tooth  variation 
in  the  isotopes  of  calculus  for  a  single  individual,  up  to  1.8%o  for  813C 
and  2.1  %o  for  815N.  In  cases  where  Poulson  et  al.  (2013)  took  mul¬ 
tiple  calculus  samples  from  a  single  individual,  they  found  a  similar 
magnitude  of  isotopic  variation.  Moreover,  this  effect  has  been 
shown  in  isotopic  studies  of  bone  when  multiple  bone  samples 
from  the  same  individual  are  analyzed  (Jorkov  et  al.,  2009;  Lucy 
et  al.,  2009),  and  for  tissue  samples  from  animals  fed  a  monoto¬ 
nous  diet  (DeNiro  and  Schoeninger,  1983).  All  of  this  suggests  that 
the  particular  bone  sampled  for  collagen  and  apatite,  and  the  exact 
tooth  chosen  for  sampling  calculus,  may  contribute  to  some  of  the 
discordance  between  bone  and  calculus  isotopic  signatures. 

The  source  of  inter-tooth  isotopic  variation  may  relate  to  the 
types  of  foods  masticated  with  molars  and  premolars  vs.  incisors 
and  canines,  and  hence,  the  organic  material  that  accumulates  in 
calculus  on  those  teeth.  It  is  not  clear  from  our  data  that  there  are 
systematic  differences  between  the  location  of  a  tooth  in  the  oral 
cavity  and  isotopic  values.  Thus,  while  three  of  four  individuals 
from  Sai  Island  showed  lower  813C  in  anterior  teeth  (incisors/ca¬ 
nines),  one  showed  elevated  values.  In  addition,  it  is  possible  that 
the  inter-tooth  differences  relate  to  the  length  of  time  represented 
by  calculus  deposits.  Thus,  calculus  on  incisors  and  canines  may 
break  off  more  readily  than  calculus  on  molars  and  premolars, 
particularly  if  an  individual  is  using  their  anterior  teeth  as  tools.  If 
an  individual  shifted  the  isotopic  composition  of  their  diet  over  the 
course  of  seasons  or  years,  this  may  be  reflected  in  the  isotopic 
composition  of  calculus  on  different  teeth.  In  this  respect,  molars 
and  premolars  may  be  more  similar  to  one  another,  and  may  be 
different  from  incisors  and  canines,  because  the  calculus  on  those 
teeth  accumulated  during  an  overlapping  window  of  time.  In  fact, 
SEM  and  TEM  studies  show  that  calculus  forms  accretionally, 
marked  by  periods  of  deposition  20-400  thick,  followed  by 
periods  of  quiescence  (White,  1997:514).  This  means  that  bulk 
samples  of  calculus  likely  represent  an  aggregate  of  different 
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depositional  episodes  and  thus  may  contain  different  dietary 
timeframes. 

Fourth,  we  suggest  that  C/N  ratios  may  be  useful  for  identifying 
calculus  samples  that  correlate  poorly  with  bone  collagen  values, 
especially  for  813C.  Empirically,  we  find  that  calculus-bone  differ¬ 
ences  are  much  less  variable  when  the  C/N  ratio  from  the  calculus 
sample  is  less  than  12.  By  contrast,  calculus  samples  with  C/N  above 
12  have  513C  values  that  are  often  highly  divergent  from  bone 
collagen,  with  calculus  813C  values  typically  highly  elevated  (less 
negative)  by  up  to  7%o.  The  effect  is  less  pronounced  for  815N.  Data 
published  in  Salazar-Garcia  et  al.  (2014)  show  a  similar  effect.  In 
their  study,  calculus  samples  with  C/N  ratios  above  12  are  also 
highly  variable  in  813C,  and  are  typically  elevated  (by  up  to  10%o). 

It  is  unclear  if  this  pattern  among  calculus  samples  with  high  C/ 
N  is  a  result  of  poor  preservation  (diagenesis),  or  represents 
reasonably  preserved  calculus  samples  that  have  a  different  bio- 
molecular  makeup,  and  hence,  unusual  isotopic  composition  rela¬ 
tive  to  bone  collagen.  We  suspect  the  former,  but  additional 
research  will  be  necessary  to  address  this  issue.  In  any  case, 
empirical  evidence  suggests  that  if  we  are  trying  to  estimate 
paleodiet  using  calculus,  and  are  seeking  to  replicate  values  that 
would  be  obtained  from  bone  collagen,  that  we  should  eliminate 
calculus  samples  with  C/N  ratios  greater  than  12. 

6.  Conclusions 

Previous  studies  on  the  C  and  N  stable  isotope  composition  of 
dental  calculus  have  resulted  in  mixed  opinions  on  the  utility  of  this 
material  for  reconstructing  paleodiet,  with  some  positive  (Scott  and 
Poulson,  2012;  Poulson  et  al.,  2013)  and  others  negative  (Salazar- 
Garcia  et  al.,  2014).  Our  independent  analyses  on  paired  calculus- 
bone  samples  from  western  North  America  and  northeastern  Af¬ 
rica  suggest  generally  positive  results.  Isotopic  signals  in  dental 
calculus  effectively  reproduce  regional-level  differences  observed 
in  bone.  These  patterns  hold  at  the  inter-individual  level  within 
sites  as  well,  but  the  correlations  are  weaker,  especially  with 
respect  to  513C  values  in  bone  apatite. 

We  suggest  that  the  poor  results  obtained  by  Salazar-Garcia 
et  al.  (2014)  relate,  in  part,  to  three  factors.  First,  they  only 
analyzed  samples  from  a  single  site  and  time  period,  the  Medieval 
Spanish  site  of  El  Raval.  We  showed  that  correlations  between 
calculus  and  bone  isotopes  are  weaker  within  sites  (vs.  across  sites) 
due  to  more  similar  diets  among  individuals,  and  hence,  the  more 
limited  range  of  isotopic  values  observed.  Second,  and  related, 
there  is  significant  inter-tooth  variation  in  dental  calculus  isotopes 
within  a  single  individual,  in  our  study  up  to  1.8%o  for  S13C  and  up  to 
2.1  %o  for  815N.  As  a  result,  the  particular  tooth  sampled  can  have  an 
impact  on  calculus-bone  isotope  correlation  coefficients,  especially 
within  a  single  site  with  a  smaller  range  of  values  for  813C  and  815N. 
Third,  calculus  samples  with  C/N  ratios  greater  than  12  are  often 
isotopic  outliers.  Including  these  samples  decreases  correlations 
between  paired  calculus  and  bone  isotope  samples  considerably. 
Removing  these  samples  from  the  El  Ravel  study  produces  results 
that  generally  mimic  those  from  bone  collagen  in  the  same 
individuals. 

That  a  positive  correlation  in  the  isotopes  of  bone  and  calculus 
exists  is  not  surprising.  Like  bone,  the  synthesis  of  calculus  begins 
with  much  of  the  same  basic  building  materials,  that  is,  substances 
that  are  put  into  the  mouths  of  individuals.  Most  of  this,  of  course, 
will  be  food  consumed  by  individuals,  though  calculus  may  be 
partially  formed  from  other  non-food  materials  present  in  the 
mouth  (e.g.,  quids;  Buckley  et  al.,  2014).  Bacteria  in  the  oral  cavity 
may  fractionate  C  and  N  in  food  differently  than  bacteria  in  the  gut, 
and  human  digestive  systems,  do.  In  other  words,  despite  the  same 
starting  substances,  C  and  N  may  be  partitioned  and  fractionated  in 


different  manners  prior  to  incorporation  into  the  structure  of  cal¬ 
culus  and  bone.  This  may  partly  account  for  the  positive  correla¬ 
tions,  but  also  the  offsets,  we  see  for  isotopes  in  dental  calculus 
(lower  813C  and  higher  815N  values),  and  may  also  add  “noise”  to  the 
bone-calculus  correlation. 

The  exact  nature  of  the  correlation  and  offset  between  calculus 
and  bone  513C  and  815N  values  is  not  well  established.  The  results 
reported  here  suggest  that  the  offset  may  vary  between  sites  and 
regions,  or  may  be  dependent  on  calculus  (and/or  bone)  preser¬ 
vation.  For  example,  calculus  813C  values  are  elevated  relative  to 
bone  collagen  at  CCO-548,  but  are  otherwise  lower  at  sites  in  this 
study.  Whether  this  variation  is  due  to  differences  in  ancient  diet, 
differences  in  oral  hygiene  practices,  the  structure  of  calculus  or 
bone  collagen,  or  preservation  conditions,  is  not  yet  clear.  In  this 
respect,  it  may  be  significant  that  CCO-548  is  the  oldest  of  the  sites 
included  in  the  study. 

What  seems  clear  is  that  further  research  is  needed  on  the 
formation  processes  of  dental  calculus.  In  particular,  we  need  a 
better  understanding  of  how  different  macronutrients  in  foods 
contribute  to  the  overall  isotopic  composition  of  calculus,  much  in 
the  way  that  we  know  that  approximately  75%  of  the  isotopic 
composition  of  bone  collagen  relates  to  dietary  protein  (Fernandes 
et  al.,  2012;  Froehle  et  al.,  2010,  2012;  Kellner  and  Schoeninger, 
2007;  Schoeninger  et  al.,  1983).  It  is  clear  that  plant  macrofossils 
and  other  organic  compounds  are  readily  preserved  in  dental  cal¬ 
culus  (Buckley  et  al.,  2014;  Hardy  et  al.,  2009,  2012;  Piperno  and 
Dillehay,  2008).  Because  plant  products  are  high  in  carbon  but 
low  in  nitrogen,  it  is  possible  that  a  significant  component  of  the 
carbon  in  calculus  derives  from  plant  products.  Because  C3  plants 
comprise  much  of  the  diets  in  California  and  Nubia,  this  may 
explain  the  lower  813C  values  we  observed  in  our  calculus  samples 
relative  to  bone  collagen.  Future  studies  could  investigate  this  by 
splitting  calculus  samples  and  treating  portions  with  particular 
chemical  washes  (e.g.,  with  acid  to  demineralize  samples,  or  with 
solvents  to  remove  organics),  focusing  isotopic  analyses  on  partic¬ 
ular  organic  fractions  of  calculus.  The  main  drawback  of  such  an 
approach  is  that  calculus  samples  are  often  quite  small  to  begin 
with  and  do  not  always  produce  enough  nitrogen  to  make  precise 
isotope  ratio  analysis  possible.  Advances  in  mass  spectrometry 
instrumentation  are  rapidly  lowering  minimum  sample  size  re¬ 
quirements  and  may  promote  such  approaches  in  the  future. 

In  any  case,  the  greater  range  of  variation  in  calculus  513C  and 
815N  values  and  the  lower  correlations  between  calculus  and  bone 
isotope  values  within  sites  is  in  many  regards  desirable.  Assuming 
calculus  isotopes  are  reflective  of  some  aspect  of  paleodiet,  it  opens 
the  door  for  complementing  rather  than  duplicating  bone  isotopic 
analyses,  which  reflect  primarily  protein  intake  in  the  case  of  bone 
collagen  and  whole  diet  in  the  case  of  apatite.  In  this  respect,  iso¬ 
topic  analysis  of  calculus  may  hold  the  potential  to  inform  on  other 
aspects  of  diet,  such  as  carbohydrate  intake.  Even  if  calculus  and 
bone  isotopes  inform  on  largely  the  same  aspects  of  paleodiet,  they 
may  still  inform  on  dietary  behaviors  over  different  windows  of 
time.  Supplementing  isotopic  analyses  of  bone  with  those  of  cal¬ 
culus  could  then  generate  additional  life  history  information  for 
particular  individuals,  and  dietary  changes  within  their  lives.  In 
short,  we  believe  there  is  still  much  potential  in  dental  calculus  to 
inform  on  the  dietary  histories  of  individuals. 
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